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Abstract 

Increasingly accurate observations of the cosmic microwave background 
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^^ ■ and the large scale distribution of galaxies necessitate the study of nonlinear 

^^ ' perturbations of Friedmann-Lemaitre cosmologies, whose equations are noto- 

riously complicated. In this paper we present a new derivation of the governing 
equations for second order perturbations within the framework of the metric- 
^ , based approach that is minimal, as regards amount of calculation and length of 

j^ I expressions, and flexible, as regards choice of gauge and stress-energy tensor. 

Because of their generality and the simplicity of their structure our equations 
provide a convenient starting point for determining the behaviour of nonlinear 
perturbations of FL cosmologies with any given stress-energy content, using 
either the Poisson gauge or the uniform curvature gauge. 

PACS numbers: 04.20.-q, 98.80.-k, 98.80.Bp, 98.80. Jk 



1 Introduction 

Over the past ten years there have been a number of significant developments in 
the theory of nonlinear cosmological perturbations and its applications that have 



* Electronic address : cl ae s . uggl a@kau.se 
T^Electronic address: jwainwri@uwaterloo.ca 



1 INTRODUCTION 2 

transformed the area from a quiet backwater into one that is at the forefront of 
research in cosmology. The main impetus for this resurgence has been the availabihty 
of increasingly accurate observations of the cosmic microwave background (CMB) 
and of the large scale distributions of galaxies. These observations necessitate the 
study of possible deviations from linearity, for example, non-Gaussianity in the CMB 
anisotropics, using nonlinear perturbations of FL cosmologie^j (see for example, 
Bartolo et al (2004a), Bartolo et al (2010b), and Pitrou et al (2010)). 

In order to motivate our work we give a brief overview of recent developments. 
First, Bartolo and collaborators have written a series of papers that apply second 
order perturbation theory to physical problems relating to the early universe and to 
the CMB. Bartolo and collaborators restrict their considerations ah initio to purely 
scalar metric perturbations at linear order and to a flat FL background. Their main 
theoretical tool is the set of expressions for the perturbed Einstein tensor at second 
order for this class of perturbed metrics, derived without restricting the gauge. 
They then introduce the Poisson gauge via gauge-fixing. We refer to Acquaviva et 
al (2003) (equation (4) for the perturbed metric and Appendix A5 for the perturbed 
Einstein tensorcl). As regards physical applications they consider the case of a scalar 
field (Acquaviva et al (2003)), a perfect fluid with linear equation of state (Bartolo 
et al (2004b)), and dust and a cosmological constant (Bartolo et al (2010a)). They 
use the expressions for the perturbed Einstein tensor referred to above to obtain the 
governing equations for the second order perturbations in these cases. 

Second, Malik and collaborators (see, Malik (2007), Malik et al (2008) and Hus- 
ton and Malik (2011)) have written a series of papers that apply second order pertur- 
bation theory to FL cosmologies with one or more scalar fields. They make the same 
simplifying assumptions as Bartolo and collaborators, but instead use the uniform 
curvature gauge. They formulate the governing equations in a way that is suit- 
able for numerical computation. Third, Nakamura (2003) introduced a geometrical 
method for constructing gauge invariants for linear and nonlinear (second order) 
perturbations which he later applied to derive the governing equations (see Naka- 
mura (2006) and Nakamura (2007)) effectively using the Poisson gauge. Finally, Noh 
and Hwang (2004) have given a comprehensive treatment of second order perturba- 
tions, with arbitrary spatial curvature and arbitrary stress-energy tensor, including 
a scalar field and a perfect fluid as special cases. They use the 3 + 1-formulation of 
the Einstein equations and write the governing equations in a so-called gauge-ready 
formO 

Despite the impressive progress that has been made to date the theory of nonlin- 
ear cosmological perturbations nevertheless presents challenges. Indeed, the theory 
is notorious for tedious calculations leading to lengthy quadratic expressions, the so- 
called source terms, that can obscure the overall structure of the equations. There 
is thus a need for a formulation of the governing equations that is both general and 
concise and that will hence provide a suitable starting point for future investigations. 



^We follow the nomenclature of Wainwright and Ellis (1997): a Friedmann-Lemaitre (FL) cos- 
mology is a Robertson- Walker (RW) geometry that satisfies Einstein's field equations. 

^See also Bartolo et al (2004a), equations (104) and (A.36-(A.43). 

•^We refer to Hwang and Noh (2007) for further details and to Hwang et al (2012) for an 
application of the formalism to determining second order perturbations of dust cosmologies. 
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Motivated by this state of affairs we present a new approach to the derivation 
of the governing equations for second order perturbations within the framework of 
the metric-based approach^ We have designed our approach to be minimal in the 
sense that we calculate the least number of objects and keep the number of lengthy 
expressions to a minimum, thereby revealing useful mathematical structure. We 
rely extensively on the notation and formalism for describing linear perturbations 
introduced in two recent papers (Uggla and Wainwright (2011) and (2012), hereafter 
referred to as UWl and UW2, respectively). In UWl we gave the linearized Einstein 
equations with an arbitrary stress-energy tensor in two different but complementary 
gauge-invariant forms, which we referred to as the Poisson form,, associated with 
the work of Bardeen (1980), and the uniform curvature form, associated with the 
work of Kodama and Sasaki (1984). In the present paper we derive and cast the 
perturbed Einstein equations at second order into a convenient form analogous to 
the linearized equations in UWl but differing by the addition of source terms. 

We begin our derivation by introducing so-called geometric perturbation opera- 
tors, which provide the first level of structure by decomposing the perturbed Rie- 
mann and Einstein tensors into a linear leading order term and a quadratic source 
term, the latter being present only at second order. A second level of structure is 
provided by the use of certain linear combinations of the perturbed Einstein tensor 
components and their derivatives, leading directly to a simple minimal set of gov- 
erning equations that is convenient for analysis. The use of a shorthand notation for 
differential operators that occur frequently and reveal key mathematical structure 
makes the equations much more tractable. A third level of structure is provided 
by our strategy of decomposing the general source terms into simpler pieces and 
identifying common expressions, without expanding them fully. Special cases of the 
source terms that appear in the literature can easily be extracted in a convenient 
form from our general expressions. 

The outline of the paper is as follows. In Section [2] we present 'geometric per- 
turbation operators' for the Riemann and Einstein tensors up to second order, and 
specialize the key operators to the Poisson gauge and the uniform curvature gauge. 
Section [3] gives the perturbed Einstein field equations to second order in both the 
Poisson and uniform curvature forms. We conclude the main part of the paper with 
a discussion in Section |H Finally Appendix |A] contains the definitions and equations 
needed to provide the background for the results in the main part of the paper, as 
well as the so-called Replacement Principle to second order, while Appendix [B] con- 
tains the detailed expressions for the Einstein source terms assuming a purely scalar 
perturbation at linear order, using the uniform curvature gauge. 



*By this we mean the standard approach to cosmological perturbations in which one formulates 
the governing equations in terms of gauge-invariant variables associated with the perturbed metric 
tensor and the perturbed stress-energy tensor, using local coordinates. 
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2 Geometric perturbation operators 

2.1 Background 

Following standard cosmological perturbation theory, we consider a l-parameter 
family of spacetimes gab{^), where gab{0), the unperturbed metric, is a RW metric, 
and e is referred to as the perturbation parameter^ We assign physical dimension 
length to the scale factor a of the RW metric and {lengthy to gab{^)- Then the 
conformal transformation 

gabie) = a^gab{,e), (1) 

yields a dimensionless metric (jabi^)- 

The Riemann tensor associated with the metric gab{^) is a function of e, denoted 
R"'^cd{^), as is the Einstein tensor, G%{e). The stress-energy tensor of the matter 
distribution is also assumed to be a function of e, denoted T%{e). We include all 
these possibilities by considering a l-parameter family of tensor fields ^(e), which 
we assume can be expanded in powers of e, i.e. as a Taylor series: 

A(e) = (o)A + e (^^A + ie^ (2)A + . . . . (2a) 

The coefficients are given by 



(°)A = A(0), «A=^ 

oe 



(2)A 



de^ 



(2b) 



where ^^^A is called the unperturbed value, '■^^A is called the first order (linear) per- 
turbation and ^^^A is called the second order perturbation of A(e). 

In particular we assume that we can expand the conformal metric gab{^) in dl]) 
in powers of e, 

gabie) = ^%ab + e (^^6 + 1^2 (2)^^, + . . . . (3a) 

We label the coefficients as 



lab:=^%ab=gab{Ql ^'^ f ab := ^'^g ab ^^'^^ 



(2)/- ._ (2)- _ 9 9ab 
J ab • yab o 2 



(3b) 

which is consistent with ( ]2bp . We refer to '■^''/afe as the /irsi order metric perturbation 
and ^^■'/afe as the second order metric perturbation. To simplify the notation, we will 
denote ^^^fab by fab when there is no risk of confusion. The conformal background 
metric 'jab and its inverse 7"'' will play an important role in that they are used to 
lower and raise indices on perturbed objects. 

2.2 Perturbation operators for the Riemann and Einstein 
tensors 

We next expand the Riemann and Einstein tensors in the form 

R^U^) = (°)i?'^^,, + e^'^R^',, + \e^ ^^^R'^\d + • • • , (4a) 

G%{e) = (°te"5 + e^^^G% + \e^ ^^^G% + ... . (4b) 



^We use Latin letters a,b, ..., f to denote abstract spacetime indices. 
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To describe the perturbations of the spacetime geometry up to second order we 
introduce dimensionless leading order linear operators and quadratic source term 
operators for the Riemann tensor R°'''cd and the Einstein tensor 

tj 6 — -K be — 2" b-ti cd- [O) 

We refer to these operators as the geometric perturbation operators and use the 
notation R"*cd(») and G%{») for the leading order operators, and 7l°'^cd{*, •) and 
Q%(^»,») for the source term operators. These operators determine the dependence 
of the hnear and quadratic terms in the Taylor series (jl]) on the perturbations of the 
metric, through equations of the following formjj 

a^(i)i?'^^ed = R^W^'^f), a'^'^R'^'a = R^U^'^f) + U'^U'^'^f, (i)/), (6a) 

a'^'^Gl = G"fe((i)/), a'^'^G% = G%{^'^f) + QM^'^f, (^V). (6b) 

It is important that there is only one leading order operator for each tensor: the same 
operator acts on both ^^'f and ^^^/. In Appendix IA.2I we derive concise expressions 
for R''''cd{^''^f) ,r = 1,2, and TZ^^di^-^^f^^^f) (see equations dUD and ([72])). 

2.2.1 Local coordinates and differential operators 

To proceed further we need to work in a coordinate frame so that we can calculate 
time and spatial components separately. We thus introduce local coordinateq3 x^ = 
(?7, X*), with ?7 being the usual conformal time coordinate! for the RW metric Qabi^)-, 
and such that the unperturbed conformal metric 7a6 := QabiS^) has components 

700 = -1, 7oi = 0, 7ij, (7) 

where 7jj is the metric of a spatial geometry of constant curvature. The function a = 
a{ri) is the background cosmological scale-factor, which determines the dimensionless 
background Hubble scalar T-L according to 

n = - = aH, (8) 

a 

where H is the true background Hubble scalar. Here and elsewhere in this paper ' 
denotes the derivative with respect to 77 of a background function that depends only 
on 7]. As in UWl we will use the geometric background scalars Ac and Cq that can 
be defined in terms of "H by the following equations: 

AG:=2{-n' + n^ + K), A'g = -{I + 3Cl)nAG (9) 

(see UWl, equation (42)). 



®Here and elsewhere we use the shorthand notation *^''-'/ for ^"^^fab, r — 1,2. 

^We use Greek letters to denote spacetime coordinate indices on the few occasions that they 
occur, and we use Latin letters i,j, fc,rn to denote spatial coordinate indices, which are lowered 
and raised using 7^ and its inverse 7'-', respectively. 

^Since we assigned a to have physical dimension length, the conformal time r/ and the conformal 
spatial line-element jijdx'^dx^ are dimensionless. We choose the x^ to be dimensionless, which 
implies that the 7^^ are also dimensionless. 
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In order to formulate the perturbation equations concisely we have found it 
helpful to introduce a shorthand notation for certain differential operators that occur 
frequently and clarify the structure of the equations. The first operator determines 
the evolution of scalar perturbations when using the Poisson gauge and is defined 
by 

£ := ^2 + 3(1 + 4)-H9^ + n^B - (1 + 3C2 )ir, (10a) 

where 

^-■=^ + ^ + 3Cg, (10b) 

(see UWl, equation (56)o). This operator has the important property that it can 
be written as the product of two first order differential operators: 

£(.) = KCaCb (^) , (11a) 

where 

CA:=dr, + nB, CB:=dr, + 2n, (lib) 

(see UWl, equation (55) and UW2, equation (39)). These operators play a central 
role, in particular when using the uniform curvature gauge. 

We will make extensive use of the second order spatial differential operators 
defined by 

D^ := 7^^D,D„ D,, := D(,D,) - '^^,^B\ (12) 

where D j denotes covariant differentiation with respect to the spatial metric 7jj . We 
will also use the following shorthand notation 

{■DAf := (D'=A)(D,A), (13) 

where A is a scalar field. 

2.2.2 Minimal representation of the perturbed Einstein tensor 

In deriving the governing equations for linear and second order perturbations in 
section [3] we follow the approach of UWl and consider four linear combinations of 
the components of the perturbed Einstein tensor and their derivatives, which we 
denote by 

(('-)(?,„ (^^te, ('-te,, ('^te°), r = l,2, (14) 

wherq^ 



(^)a;. := TM.^'^G'^.) - h^/'^G^, (15a) 

4(''te°o + l^"^G\, (15b) 

-Di^'te^o - 3H("te° , (15c) 



Hj 



^"^Gi 



^See the final paragraph of section 3.2 in UWl for references to the literature where related 
forms of this operator appear. 

-"^"Note that ^'^^Gij = ^kS'^'^G^j is not symmetric when r ~ 2. We thus compensate by symmetriz- 



ing in (|15ap . so that ^^''Gij, as defined by that equation, is symmetric. 
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and Cq is defined by ([9]). The use of these combinations leads to a convenient form 
of the governing equations. In particular, in the linear case (r = 1) when using the 
Poisson gauge the expression fll5bp leads directly to the evolution equation for the 
Bardeen potential in terms of the differential operator C in flTO|) . while fll5cp leads to 
the generalized Poisson equation that determines the density perturbation in terms 
of the Bardeen potential (see UWl equations (54)). We will see that this pattern is 
repeated in the nonlinear case. 

In order to calculate the Einstein combinations (TT5|) we express them in terms 
of the perturbations of the Riemann tensor using ([5]) and 



3xpi 





^'^G,,=lH^''^R%\j), (16a) 

('■te = -|(1 + ^Clt^R%, - l^'^R%k, (16b) 

('■te, = iD,M/2P% - Sn^'^R'^^k, (16c) 

It follows from fl6b|) that the dependence of the Einstein combinations (llSp on 
the metric perturbations is given by equations of the form 

a2(2)G,, = G,,p)/)+^,,(/,/), (17a) 

(^)G = Gp)/)+^(/,/), (17b) 

(2)G, = G,((2)/)+^,(/,/), (17c) 

a^(^)G''. = G° ((2)/)+^° (/,/), (17d) 

where for brevity we have denoted '■^^/ by / in the source terms. 

2.2.3 The leading term operators 

We can express the leading order operators in ( IT71) in terms of R"'^cd{f) using equa- 
tions (1161): 



a^^^V?., = 


= G..™, 


a^dte = 


= G(«/), 


a2(i)G, = 


= G.(«/), 


a2(i)G° = 


= G°(W/), 






G.,(/) - 


= 7fe(iR"1a|i)(/), 


(18a) 


G(/) = 


= -i(l + 3C^)R%,(/)-|R°V(/), 


(18b) 


G.(/) - 


= |D,R%,(/)-3HR°^,,(/), 


(18c) 


G°.(/) 


= R°'..(/), 


(18d) 



"'^"'^Here and elsewhere (..) stands for a symmetrized and trace-free spatial index pair of a perturbed 
object with respect to the background spatial metric "/ij, while objects with a symmetrized and 
trace-free spatial index pair are given a hat over the kernel symbol. 
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where / denotes ^^^f or ^'^^f. On using equations (17T|) for R'^^cd(/), in conjunction 
with (TMl) . we obtain after some manipulationo 



(19b) 



G,,(/) = iD,,(/oo - |/\) + /:bY,,(/) + D'^^,/,), - |(D2 + 3ir)4-, (19a) 

G(/) = -i [(£ - C'a^')f\ + {S-HCa + D^) (f,, - If, 
- 1(1 + 34)D^„n + f (-Cs + 3^4) DVfco, 

G,(/) = -Di [i(D2 + 3K)f\ - iD^„n] + SHD'^ Y,,(/), (19c) 

G°(/) = fD,Y(/)-D^Y,fc(/), (19d) 

where 

Y(/) := Infoo + ldj\ - DVfco, Y,,(/) := IdJ,, - D(,/,)o, (20) 

and / denotes ^^'f or ^^7. These equations are the main resuh of this subsection 
and play a central role in our derivation of the perturbation equations. 

2.2.4 The source term operators 

Calculating the Einstein source terms Q{f, f) in (TT7I) presents a major challenge. 
Our strategy is to first express them in terms of 'R-"'^cd{f, f) using equations (fT6!) : 

QMf)=7kiin'^%\j)ifJ), (21a) 

g{f, f) = -|(i + 3C2 )7^%,(/, /) - |7^°V(/, /), (21b) 

G,{f, f) = iD,7^%,(/, /) - snn'^kif, /), (21c) 

Q%{fJ) ='R'\k{fJ). (2id) 

We then substitute the expression for 7l"'\d{f, f) given by equations fl72l) in the 
Appendix into fl2Tl) . This leads to the expressions f l75|) for the Einstein source terms, 
each as a sum of simpler terms that can be calculated separately. The constituent 
terms are given by equations (1761) -( 178|) . 

2.3 Gauge invariants for the Einstein tensor 

We associate gauge invariants with the linear and second order perturbations of 
any tensor using a method pioneered by Nakamurao, which we modify to ensure 
that the gauge invariants are dimensionless. The process of construction, which 
involves introducing co-called compensating gauge fields denoted by ^^^X, r = 1,2, 
is described briefiy in Appendix lA.Sj and is specified by equation (18 ip for an arbitrary 
tensor and by f E7|) for the conformal metric tensor. The gauge invariants associated 
with the perturbations of the conformal metric tensor and the Einstein tensor by 
this process, which we refer to as X-compensation, are denoted by *^^''fab[X] and 
^^) G% [X] , respectively. 

In equation (I6bl) we expressed the perturbations ^^"^0%, r = 1,2, of the Einstein 
tensor in terms of the perturbations ^^'fab of the metric tensor in gauge- variant form. 



-^^We refer to footnote [TT] for the (..) and hat notation, and to equations pll |) - p^ for the 
definitions of the differential operators. 

^■^See for example Nakamura (2007), section 2.3. 
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using the Einstein geometric operators G^ and Q%. By applying the Replacement 
Principle^ 1 to equation (]6bp we obtain the following expressions for the Einstein 
gauge invariants in terms of the metric gauge invariants: 

(i^C^^X] = G'^,((i)f), (2)Qa^[X] = G'^b((2)f) + g%{WfWf), (22) 

where we use the shorthand notation "^^^f = '^'')faft[X], r = 1,2. 

In this section we derive explicit expressions for the gauge invariants on the 
right side of (122]) for two specific choices of the gauge fields X. The starting point 
is to consider the gauge invariants associated with the first and second order metric 
perturbations. 

2.3.1 Gauge invariants for the metric tensor perturbations 

To construct dimensionless gauge invariants associated with the linear perturbation 
fab of the conformal metric tensor we definq^ 

iab[X] := fab - a~'^£wx (a^7ab) , (23) 

where the gauge field (^)X" has to be chosen appropriately. In order to construct a 
metric gauge field one has to decompose fab into scalar, vector and tensor modes, 
which we label as follows (see UWl, equation (18)): 

/oo = -2^, (24a) 

foi = BiB + Bi, (24b) 

fij = -2^7,, + 2D,D,-C + 2D(iC,-) + 2Q,, (24c) 

where the vectors Bi and Ci and the tensor Cij satisfy: 

n'Bi = 0, WCi = 0, C\ = 0, Ti'Cij = 0. (24d) 

We use equations fl2^ as a model for doing a mode decomposition of fa6[X], using 
an obvious notationO 

As shown in UWl there are two ways to choose X uniquely in terms of fab, 
leading to the Poisson gauge field Xp and the uniform curvature gauge field Xc. 
The corresponding expressions for f^fefX] are as follows (see UWl, equations (28)- 
(31)): 

f,,[Xp] = -2^7,, + 2Q,-, (25a) 

* := nX,], (25b) 

BiB + Bi, f,,-[Xe]=2Q„ (26a) 

i^Make the replacements (')/ -> (''^f [X], r := 1, 2, in ^, in analogy with (|83l). 
^^See equation ()87a|) in this paper and equation (16) in UWl. 





foo[-^p] = —2$, foj[^p] = Bj, 


where 






$ := $[Xp], 


and 






foo[-^c] = —2 A, foj[^c, = 



16 



For example, ip — > ^[X],B -^ B[X], as in equation 
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wherq^ 

A:=^X,], B:=B[Xe]. (26b) 

In both case the vector and tensor modes satisfy 

WBi = 0, C\ = 0, WCij = 0. (27) 

To construct gauge invariants associated with the second order perturbation ^'^^fab 
we introduce a second gauge field "^^^X and definq^ 

^Xb[X] := (2)/„, - a'^£c2,j, {a^^ab) + J'abiX], (28a) 

where 

Tab[X] := -a-^£wx{'2a^fab - £wx{a^lab)). (28b) 

The key point is that one can construct a gauge field "-^^Xp such that "^^^afofXp] has the 
same form as fa6[Xp] in (!25l) . In other words, '^^¥af,[Xp] can be obtained by making 
the substitutions 

$^(2)$^ vl/^(2)vl>, B,->(2)B,[Xp], Q,^(2)Q,[Xp], (29) 

in (j25]). Similarly, one can construct a gauge field '^^-'Xc such that *^^^ab[Xc] is obtained 
by making the substitutions 

A ^(2) A, B^(2)B^ B,^^'^B,[X,], Q,->(2)Q,[X,], (30) 

in 0261) . Details about the construction of these gauge fields and the expressions for 
the metric gauge invariants in equations (123]) . (12^ . (12^ and O30p in terms of the 
gauge-variant metric perturbations are given in Appendix IA.4[ It is important to 
note, however, that these explicit expressions are not required in what follows. All 
that is required is the general form of fab[X] and '^^¥af,[X] in the Poisson gauge and 
in the uniform curvature gauge, as given by equations (1251) . fl2B]) . fl2Ul) and fl5U]) . 

2.3.2 The leading order terms 

To obtain gauge-invariant expressions for the leading order terms we simply make 
the substitutions fab -^ fa6[X], ^'^^fab — ^ '•^■fafefX] in equations (fT9|l . For the Poisson 
gauge we use (l25l) which gives the leading order Einstein operator G acting on the 
first order metric perturbation fp = fa6[Xp]. After some manipulation we obtain^ 

Q,(fp) = JDiji^ - $) - D(,/:bB,) + {Csdr, + 2K- D^) Q,-, (31a) 

G(fp) = 2 [(£ - 4d2) M/ - {HCa + |d2) (vl/ - $)] , (31b) 

G,(fp) = 2D,(D2 + 3K)^ - ln{T>^ + 2K)B„ (31c) 

G° (fp) = -2D,(9^^ + H^) + i(D2 + 2K)B,. (31d) 



-'^'^In UWl we introduced the symbols A and B for these gauge invariants, foUowing the notation 
of Kodama and Sasaki (1984), equations (3.4) and (3.5). 
^^See equation dHZ]) in Appendix IX4l 
i^The identities (B.39b) and (B.39f) in UWl are needed. 
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To obtain G acting on the second order metric perturbation ^^^p = ^^^^^[Xp] we 
simply make the replacements ( 129|) in ( 13T]) . 

Similarly for the uniform curvature gauge we use fl26p . which gives the leading 
order Einstein operator G acting on the first order metric perturbation fc = faf,[Xc]. 
After some manipulation we obtaino 

GijiQ = --DijiCeB + A) - Bf^iCeBj) + (£^9^ + 2K - D^) Q„ (32a) 

G(f,) = 2n[CAA + C^D'B] + |D2(£bB + A), (32b) 

Gi(fc) = -2-HDi(D2 + 3K)B - ^HCD^ + 2K)B,, (32c) 

G°i(fc) = -2Bi(nA - KB) + i(D2 + 2K)Bi. (32d) 

To obtain G acting on the second order metric perturbation ^^^c = '■^^af,[Xc] we 
simply make the replacements ( 130|) in ( 132|) . 

Equations ( 13T|) and ( |32|) provide the leading order terms in the perturbed Einstein 
equations (l40l) at linear and second order using the Poisson gauge and the uniform 
curvature gauge, respectively. For the reader's convenience we note that the various 
differential operators are defined by equations (IT0l)-(IT2l). 

2.3.3 The source terms 

The source terms are given in general in gauge- variant form by fITS]) in conjunction 
with fl751) - fl75|) . They are obtained in gauge-invariant form by simply making the 
replacement fab — ^ fabl-^] in these equations, in particular with fab[Xp] for the Pois- 
son gauge and fa6[Xc] for the uniform curvature gauge (see equations (125|) and (!26l) ). 
To illustrate our approach we consider a popular special case, namely the Poisson 
gauge with the metric perturbation restricted as follows: 

Metric assumptions: The vector and tensor modes of the metric perturbation are 
zero at first order, i.e. 

Bi = 0, Qj = 0, (33a) 

and in addition the scalar mode at first order satisfies 

$ = ^. (33b) 

Subject to these assumptions fl25l) reduces to 

foo[Xp] = -2^, ioi[Xp] = 0, fij[Xp] = -2^jij, (34) 

where "^ is the Bardeen potential. We now make the replacement fah — ^ fab[-^p] in 
the expressions fl75|) for the Einstein source terms using the special metric perturba- 
tion flMl) . The constituent terms, as given by equations fl76|) - fl78|) . can be evaluated 
separately and one finds that many terms are zero. This calculation yields the 
following simple expressions for the source terms: 

^,,(fp,fp) = 4(D,,^2 _ (D^,^)(D,)^)) , (35a) 

^(fp, fp) = -|(1 + 34) 7^(^, ^) - |(D^)=^ - mCA'^^ (35b) 

^,(fp, fp) = D,7^(^, m) + 127^(a,^)D,^, (35c) 

^° (fp, fp) = m-D,^^ - 4(a,^)(D,^), (35d) 



""The terms involving the vector and tensor modes are the same as in the Poisson case. 
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where 

7^(^, ^) := 2 [3(9^^)2 - 5(D^)2 + 4(0^ + 3^)^^] . (35e) 

Note that the source terms are quadratic expressions in the Bardeen potential \I/ and 
its derivatives. Nakamura (2007) has given the source terms in this case in the form 
Q%{ip, fp). We find complete agreement with his equations (6.13)-(6.16) when they 
are transformed into the form (j35])- 

In summary, equations (!3T|) and ( 135|) give the constituent geometrical parts in 
the second order field equations P0|) in section [31 They thus provide the foundation 
for determining second order metric perturbations in the Poisson gauge, subject to 
the simplifying metric assumption fl33p . 

Finally we note that one can similarly derive expressions for the source terms 
using the uniform curvature gauge, subject to the simplifying assumption f l33ap . 
The resulting expressions are more complicated than f l35|) and so we give them in 
Appendix [Bl 

3 Perturbed Einstein equations 

3.1 General structure of the governing equations 

At zeroth order the non-zero components of Einstein's field equations are given bjo 



a 



a^ 



(°te°o = -3(?/' + K) = -a'^% =a2(o)r°o, (36a) 

^''^G) = -{2n' + H^ + K)6) = a^'^yS) = a^^'^'T), (36b) 



where "H is given by ([H]) and K is the curvature index, defined in equation (jST]) in 
Appendix lA.ll 

The perturbed Einstein equations at linear and second order are given by 

(i)^a^ = {^)T%^ mG% = ^^^T%. (37) 

Assuming that the background Einstein equations are satisfied we can write these 
equations in terms of the gauge invariants '^^'•'G'^fX] and '^'')T'fe[X]: 

(1) G"4X] = (1) T% [X] , (2) Qo^ ^x] = (2) T", [X] , (38) 

as follows from the definition flHTj) in Appendix IA.3I We can now use f l22|) to express 
the left side of these equations in terms of the Einstein operators: 

G»,((i)f) = Wt%[X], (39a) 

G"b((2)f) + ^'^,((i)f,(i)f) = (2)T",[X]. (39b) 

At this stage we are considering an arbitrary stress-energy tensor, whose components 
are regarded as primary objects, i.e. they are not constructed from other quantities 



^^See, for example, Mukhanov et al (1992), equation (4.2), noting the difference in signature. 
We use units c = 1 and 8nG = 1, where c is the speed of hght and G is the gravitational constant. 
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as in the case of the Einstein tensor. Before continuing we note that equations 
correspond to equations (2.49) and (2.50) in Nakamura (2007) o 

We next consider the combinations of equations fl39|) corresponding to the com- 
binations of the Einstein components defined in flTSl) : 

G,,((^)f) = (^)t,,[X], G,,((2)f) + ^..(f^ f) = (2)t,,[X], (40a) 

G((i)f) = (i)t[X], G((2)f) + ^(f, f) = (2)t[X], (40b) 

Q(Wf) = «T,[X], Q((2)f) +^(f,f) = (2)T,[X], (40c) 

GO((i)f) = (i)tO[X], G°((^)f) + e^O(f,f) = (2)tO[X]. (40d) 

The hnear combinations of '^'')T'b[X] in (1401) are defined in analogy with (TT^ b 

(-)T[X] 



7.(.('->r*^,)W-h/^>r\[X], (41a) 

c|(^>rOo[x] + |('->r\[x], (41b) 

-D,('->r°o[X]-3?^('->r°[X], (41c) 



where 



C| = ,/-, ^y=\^'^T% ^% = -('yT'o. (42) 



(o)p' 



Equations (HOj) give a convenient minimal form, of the governing equations for 
linear and second order perturbations for any choice of gauge field X and any stress- 
energy tensor. The leading order terms G(f), where f = "^^^f or *^^)f, are given 
by (PTj) for the Poisson gauge X = Xp and by (152]) for the uniform curvature gauge 
X = Xc- The source terms Q{{, f ) are obtained in general by making the substitution 
fab — )■ fab[X] in equations fl75|) - fl78|) for an arbitrary X, and are given directly by f l35|) 
for the special metric perturbation fl33|) when using the Poisson gauge. 

The final step is to decompose the governing equations PUI) into equations for 
the scalar mode, the vector mode and the tensor mode. As with the metric we 
perform a mode decomposition of the stress-energy gauge invariantso 

(^'■>r,, [X] = D, >)n + 2D(>)n,) + ^'%, (43a) 

(^^>r,[X] = D,MA + ('■)A„ (43b) 

('•>r°[X] = D,(^V + (^')V^i, (43c) 

(">r[x] = (">r, (43d) 

where 

JD'^'^li = 0, ^'\l\ = 0, Bi^'\l) = 0, B'^'^Ai = 0, B'^'K'i = 0. (44) 



^^See also, equations (38) and (39) in Nakamura (2006). Nakamura's metric gauge invariants 
are related to ours according to Cab ~ a? ^^^iab[X], Hab = c? '^'^^iab\X\. 

^"^At first order '^^T'ij[X], *^^-'T[X], '^^)Ti[X] are intrinsic gauge invariants since they do not 
depend on the choice of gauge field and hence can be written as ^^T'ij, ^^^, ^^^Ti (see UWl, 
section 2.3). At second order this is no longer the case. 

^^For brevity we drop the argument \X\ for the various mode terms. 
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A difficulty arises that is not present at the hnear level. The leading order terms 
G(*^^)fp) and the stress-energy terms *^^^T[X] in (HOj) are expressed explicitly as a sum 
of a scalar term, a vector term and a tensor term. On the other hand, the source 
terms ^(f, f) do not have this form, as can be seen, for example, from (135|1 . One 
thus has to apply what we call mode extraction operators to ( 1401) in order to separate 
the modes in the source terms. These operators are defined in (1851) . using the letters 
S, V and T to denote the scalar, vector and tensor modes, respectively. 

For later use we note that applying the mode extraction operators to (H3l) yields 

ir)n = s'^(^%[x], (^)n, = v/'^(^>r,4x], (^)n,, = T^r^'%g[xi (45a) 

(r)^ = 5i(r-)x . [X] , ('■)Ai = V/(''>rj [X] , (45b) 

{r)y ^ 5i(r-)xO. [X] , (■'^- = V^^^'^^j [X] . (45c) 

3.2 The mode-decomposed governing equations 

In this section we give the mode-decomposed form of the governing Einstein field 
equations at second order for perturbations of an FL cosmology with arbitrary mat- 
ter content, ffist using Poisson gauge invariants, and then using uniform curvature 
gauge invariants. The source terms, identified by the kernel Q, are obtained by 
making the substitution / — )• fp or / — )• fj. in equations (!75|) -( [78|) . or directly by ( 135|) 
for the special metric perturbation (!33l) when using the Poisson gauge. We note 
that, in accordance with ( HOl) . the governing equations at first order can be simply 
obtained from the equations at second order by dropping the source terms, indicated 
by the kernel Q, and dropping the exponent '•^^ (or replacing it by '•^^). 

The Poisson form 

To obtain the Poisson form we substitute (|3T|) and (H3l) into (j40|) and then apply 
the mode extraction operators (|85|) . which leads to: 

Scalar mode 
(£-C^D2)(2)^ = i(2)rp+ 

{hCa + p^) ((2)np - s'^ ^,,(fp, fp)) - |^(fp, fp), 

(D^ + 3ir)(2)vl/ = 1(2) Ap - IS' e^,(fp, fp), (46c) 

Vector mode 

CB^'^^[X^] = -2(%4Xp] + 2V/' ^,fc(fp, fp), (47a) 

(D^ + 2ir)(2>B,[Xp] = 2(2)K[Xp] - 2V/ ^°-(fp, fp), (47b) 

Tensor mode 

{CBd, + 2K-BY^C,j[X^] = (%.[Xp] -7^'="^fc„(fp,fp). (48) 



(46b) 
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Here *^^^\1/, (^)$, (^fej[Xp] and '•^^Cjj[Xp] are the second order Poisson metric gauge 
invariants and fp is shorthand for the first order metric perturbation fQfe[Xp] as given 

hym- 

The evolution of the scalar perturbations is governed by equation (146bp . a second 



order partial differential equation for (^^\E'. In order to obtain a solution one first 
has to solve the linearized field equations for the first order gauge-invariant metric 
perturbation fafe[Xp], which then determines the Einstein source terms. Once the 
second order matter terms '■^-Tp and '^^-Tlp have been specified, one can solve (I46bp 



for (^)\l/ and then successively use (146 ap . f l46cp and ( ]46dl) to calculate (^^$, ^^^Ap and 
^'^^Vp, respectively. 

The uniform curvature form 

To obtain the uniform curvature form we substitute (!32|) and ( H3|) into (HOj) and then 
apply the mode extraction operators. For the sake of brevity we give only the scalar 
mode since the vector and tensor modes have essentially the Poisson form. 

Scalar mode 

£s(2)B + (2) A = -(2)n, + S'^ g,j{f,, fj, (49a) 

n{c/'^A + clB'^'^B) = |(2)r^ + iD2[(2)n,-5^^-^,^.(f,,f,)]_i^(f„f,), (49b) 

•H(D2 + 3ir)(2)B= -i(2)Ae + |5*^,(fe,fe), (49c) 

7^(2) A _ K(2)B = -i(2)K + ^S' ^° (fe, fe). (49d) 

Here *^^^A and '■^^B are the second order uniform curvature metric gauge invariants 
and fc is shorthand for the first order (in time) metric metric perturbation faf,[Xc] 
as given by ( l26|) . 

These equations differ from the Poisson form ( H6|) in that the evolution of the 
scalar potentials A and B is governed by two coupled first order partial differential 
equations (149 ap and (]49bp . In order to obtain a solution one has to follow a two 



step procedure, as with the Poisson form. 

Commentary 

The systems of equations (l46l) and (l49l) . together with the expressions (ff5l) -(j78ll 
for the source terms, are new and constitute one of the main results of this paper. 
Either set of equations determines the behaviour of second order scalar perturbations 
of an PL cosmology with arbitrary stress-energy content. We emphasize that the 
specific form of the evolution equations P6l) and fH9|) depends on the matter terms, 
specifically, the P-terms and the H-terms, both at first order and at second order. 
These quantities are determined by the stress-energy tensor using equation fl84bp in 



conjunction with fl43dp and (145 ap . In the next section we illustrate how to calculate 
these quantities for the simple case of a perfect fluid. 

Equations (H6l) and (H9!) illustrate a fundamental difference between second order 
and first order perturbations. The analysis of linear perturbations is simplified by the 
fact that the three modes, namely scalar, vector and tensor, decouple and hence can 
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be analyzed separately. At second order each mode involves leading order terms, 
described by the same operators that occur at first order, but also complicated 
quadratic source terms that, in general, contain all linear modes, which leads to a 
phenomenon one may refer to as source mode coupling. For example, this means 
that the first order metric perturbation fp that determines the Einstein source terms 
in (1461) contains scalar, vector and tensor modes in general (see (123]) ). In other 
words the vector and tensor perturbations at linear order contribute to the scalar 
perturbation at second order. On the other hand equations fH6l) show that a purely 
scalar linear perturbation [i.e. if the vector and tensor modes at the linear level are 
assumed to be zero, as is often done) generate all three modes at second ordero 

3.3 Matter gauge invariants for a perfect fluid and A 

In this section we determine the second order stress-energy perturbations (^^[X], 
(^tl[X], *^^tlj[X] and ^^tljj[X] that appear in the governing equations in Poisson form, 
as given by equations (H6|) . fHTj) and (HHl) . when the stress-energy tensor describes a 
perfect fluid and a cosmological constant A: 

n = (p + p)u''ub + p6% - An. (50) 

Here u", p, p are the fluid's 4- velocity, energy-density and pressure, respectively. 
Since A is e-independent it follows that it does not appear in ^^^% and ^'^^%, and 
it only affects the perturbed field equations indirectly via background quantities 
determined by the zeroth order field equations. For simplicity we assume p = wp 
with w = constant, which implies that 

(^^p = C^(^^p, with C^ = w, for r = l,2. (51) 

In order to find the desired quantities we need to calculate "^^^[X] and ^^^jj[X] 
for the stress-energy tensor (15n|) . We begin by expanding ( ISUjl to second order, which 
yields 

a' 



,2 I 



(^n = -a' (^)p - 2AW^^% {^% - /o,) , (52a) 

^ (2>r\ = 3a' (^)p + 2AtY^'^\ {^% - /o.) , (52b) 

a'^'% = 2AT{^%-foii}^%, (52c) 

where ^^^Vi := a~^^^\i. Applying the Replacement Principlcj gives 

^'yr'o[X] = -%[X]-UMX], (53a) 

(2>r\[X] = -3C^(2)p[X] + UMX], (53b) 

(2>r,,[X]=W,,[X], (53c) 

where 

W,,[X] := 2At'^'V,[X] {^%[X] - ArioAX]) , (54) 



^^See, for example, Bartolo et al (2004a), section 3.1. 

^^To obtain (^'T'^bi^] replace each perturbation variable on the right side by its gauge invariant 
formed by X-compensation i.e. ^^^ — )■ ^^^p[X] and Ar^^'^Vi ~¥ ^^X^i[X]. See equation ([84)) . 
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and 

At := a\(y+^y), ^'M[X] := (i>r°[X]. (55) 

Equations fl53al) and fl53bp . together with (I41bp . result in 

(2)r[X] = |(l-3C2)W\[X]. (56) 

As expected, we see that "^^^[X] and *^^'Tjj[X] are purely source terms 1^1 

For simplicity we now assume that the linear vector modes Vi[X] and Bj[X] are 
zero, i.e., 

V,[X] = B.y[X], fo,[X]=D,B[X]. (57) 

We then choose X = Xp, noting that B[Xp] = and V^[Xp] = V. It follows that 

W,,[Xp] = 2AT\'D^V){B,V), (58) 

and hence 

(2)r[X] = 1(1 - 3C^)AT\rtVf, ^^^iAXp] = 2^^iD(iVD,)K (59) 

The required quantities (^tl[Xp], *^^tlj[Xp] and ^^tli^fXp] are then obtained by ap- 
plying the mode extraction operators as in (I45ap . 

In order to facilitate comparison with the literature we relate the gauge invariants 
associated with the matter density and with the matter velocity in the Poisson gauge 
at second order, denoted ^^-"Pp, and ^^Vp, respectively, to our stress-energy gauge 
invariants ^^■'An and ^^■'Vn: 



^% =(2)Ap + m^'\ - 

Pp = A + 3H\/p, 


2A^\'DV^)\ 

\( Pr. \ 1 


J 

) 


(60a) 
(60b) 

(60c) 





where 



3.4 The governing equations in a simple example 

We now specialize the general governing equations P6l) . (HTI) and (148|) at second 
order in the Poisson gauge to a perfect fluid and a cosmological constant with only 
scalar first order contributions and $ = \E', i.e. we impose the metric assumptions 
(155]) . We hence substitute for the metric source terms ^(fp,fp) from (1551) and for 
the stress-energy perturbations '•^^P and the three ^^^Il-terms from (145 ap and (I59j) . 
This leads to the following governing equations: 

Scalar mode 

(% _ (2)$ = _4^2 ^ 2S'^Mij{^, ^), (61a) 
{C - C^D2)(% = 1(1 - ?>Cl)A-^\T^Vf + 1(1 + 34) 7^(M/, ^) 

+ |(D^)2-|D2^2^2(-H£^ + iD2)5'^.M,,(^,^), ^ ' 

(D^ + 3K)(% = i(2)A - i7^(^, ^) - QHS\{dr,^){p,^)) , (61c) 

{dr, + H)(% = -i^'V - 2S\{dr,^){Tii^<)) + 2nS''Mij{^, ^), (61d) 



^^Since (^^[X] — and '^^^^^[X] = for a perfect fluid it folfows that the feading order term in 
(2)r[X] and in ^^^.,j[X] is zero. 
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where 

n{-^, ^) := 2 [3{drj-^f - 5(D^)=^ + 4(D=^ + 3K)^'^] , (61e) 

M,,{^, *) := 2(D(,^)(D,)^) + A^\T>(,V){B^)V), (61f) 

and V is given bjo 

^ = -2(5, + 7/)^. (61g) 

Vector mode 

C^i^JB, = -4V/^A^,,(^, ^), (62a) 

(D^ + 2i^)(2)Bi = 2(2V. + 8V/((9^^)(D,-^)). (62b) 

Tensor mode 

{Csdr, + 2K- D^) (2)Q^. = 27-/'"A^fc^(*, ^). (63) 

Commentary 

An attractive feature of the Poisson gauge at first order is that if the anisotropic 
stress is zero (for example, for a perfect fluid or a scalar field) then $ = \E', i.e. 
there is only one metric gauge invariant for the scalar mode. Equation (I61al) shows 
that this feature is not preserved at second order due to the presence of the source 
terms. There are thus two second order metric gauge invariants, '•^^$ ^ ^^)\&. We 
shall refer to *^^-'$ as the second order Bardeen potential and to ^^^^ as the second 
order Bardeen curvature^ to distinguish their roleso 

The system of equations (l6Tl) - (!63l) is closely related to but simpler than equations 
that appear in the literature. Nakamura (2007) has given a system of equations 
that can be transformed intqfj our equations (!6T|) - (!63|) . There are, however two 
important differences. First, Nakamura chooses the Bardeen potential *^^^$ rather 
than the Bardeen curvature ^^Hl/ to be the primary metric gauge invariant for the 
scalar mode, which has a major drawback: the source terms contain second order 
time derivatives and are significantly more complicated. The second difference is in 
the treatment of the density perturbation. We use '^^•'Ap, which satisfies a generalized 
Poisson equation and can be viewed as being analogous at second order to the well 
known Bardeen gauge invariant e^- In contrast, Nakamura uses the gauge invariant 
^^•"Pp associated with the density perturbation at second order in the Poisson gauge, 
as given by ([60a]) cj 

Equations (!62|) and (!63l) show that a purely scalar linear perturbation gives rise 
to a vector and a tensor perturbation at second order, with the link provided by 
the tensor Aiij. The gravitational waves described by this tensor perturbation have 
been investigated in detail by Ananda et al (2007) and Baumann et al (2007). 



^®This is obtained from the governing equations for the first order perturbations. 

^^Strictly speaking one should also make this distinction at first order, but since $ = ^ in many 
applications we simply refer to \1/ as the Bardeen potential. 

■^"His equations (6.38), (6.44), (6.4f) and (6.42) can be used to obtain the first four of our 
equations (|6ip . his equations (6.39) and (6.33) can be transformed into our equations (|62|) . and 
finally his equation (6.40) yields our equation (I55|) . We also refer to Nakamura (2006) for a brief 
summary. 

^-^Note that ^^^p[Xp] = o? ^^^e in Nakamura's notation. 
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4 Discussion 

The systems of equations (1461) and (H9l) which govern nonhnear perturbations, to- 
gether with the expressions (175l)-( I78l) for the source terms, are new and constitute 
one of the main results of this paper. Because of their generahty these equations 
provide a starting point for determining the behaviour of nonhnear perturbations of 
FL cosmologies with any given stress-energy content, using either the Poisson gauge 
or the uniform curvature gaugeo These equations exhibit the same concise struc- 
ture as the governing equations for linear perturbations given in UWl (see equations 
(52) and (54)), in which the evolution of the metric perturbations is determined by 
the second order factored differential operator C in the Poisson form and by the 
pair of first order differential operators Ca and £b in the uniform curvature form. 
This structure arises directly from our use of specific linear combinations of the 
components of the perturbed Einstein tensor and their derivatives, as given by (llSp 
(see UWl, equations (39) and (40) for the motivation). Indeed the three operators 
are visible at an early stage in the derivation in the leading order expressions flT9|) 
for the perturbed Einstein tensor. This is in contrast to the literature where it is 
customary to simply calculate all the components ^'^'^G% of the perturbed Einstein 
tensor, and then form linear combinations of the perturbed Einstein equations (see 
for example, Acquaviva et al (2003), equation (4) for the perturbed metric and Ap- 
pendix A5 for the perturbed Einstein tensor, and Bartolo et al (2004a), equations 
(104) and (A.36)-(A.43).). This process involves more extensive calculations than 
in our approach and may lead to expressions that are not optimally simplified while 
hiding important mathematical structures. 

There is one issue that deserves particular attention, namely the fact that there is 
no unique choice of gauge invariant associated with the perturbations of the matter 
density. At the linear level there are three commonly used choices, the Poisson 
gauge invariant, the uniform curvature gauge invariant that is related to one of the 
so-called conserved quantitieo, and the total matter (or co moving) gauge invariant. 
The last-mentioned is the well-known Bardeen gauge invariant, and is related to the 
spatial gradient of the matter density orthogonal to the fluid flow. At the second 
order level the situation is more complicated and requires further invest igat ion o 
The appropriate choice may depend on the physical situation under consideration. 

In this paper we have focussed exclusively on using the perturbed Einstein field 
equations to describe the dynamics of nonlinear perturbations. There are, however, 
two alternatives to the direct use of the Einstein equations. First, one can use the 
perturbed conservation equations for the stress-energy tensorjf^ and second, one 



^^ For example, by specializing these equations we can derive in an efficient manner the various 
equations for second order perturbations that appear in the papers by Bartolo and collabora- 
tors (Poisson form) and by Malik and collaborators (uniform curvature form), referred to in the 
introduction. 

^^See, for example, Malik and Wands (2004), equations (4.17) and (4.27). 

^"^See Bartolo et al (2010a), section 3 and in particular equation (29) for '-^W in a AC DM model, 
Christopherson and Malik (2009), equation (4.8) for ^^^p^ and Noh and Huang (2004), equation 
(273) for (2)p^. 

■^^See, for example Bartolo et al (2004b), equation (4.3), and Noh and Hwang (2004), equations 
(104) and (200). 
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can use the 1 + 3 formalismjfj expanding the exact equations to second order and 
making them gauge-invariant. More work needs to be done in this regard. An 
additional aspect of the dynamics of scalar perturbations that we have likewise not 
touched on is that under certain conditions [i.e. in the long wavelength regime) 
the governing equations admit so-called conserved quantities, i.e. quantities that 
remain approximately constant during a restricted epoch. These quantities, which 
were initially introduced for linear perturbations (see, for example UW2 section 4 
for a unified overview) have now been generalized to second order perturbations J^l 

Acknowledgments 

CU is supported by the Swedish Research Council (VR grant 621-2009-4163). CU 
also thanks the Department of Applied Mathematics at the University of Waterloo 
for kind hospitality. JW acknowledges financial support from the University of 
Waterloo. 



A Derivation of the perturbation equations 

A.l Exact curvature expressions 

In UWl we derived an exact expression for the Riemann curvature tensor R"-\ci{e) 
of 5'afe(e) = o^Qabi,^) by replacing the covariant derivative "Va of Qabi,^) with the 
covariant derivative °Va of QabiS^) = lab- We first make the following definitions: 

ra:=°Va (In a), (64a) 

QUe) := r%) (%^c)d(e) - |°Vrf^,,(e)) . (64b) 

Since ^ajbc = it follows that 

Q%M = 0. (65) 

The desired expression is as follows (see UWl, (B.8), (B.IO) and (B.12b)): 

a^R'^Ue) = R^'cdie) + A5\u\{e), (66a) 

where 

R'^'cdie) = -g'\e) (^R\^ + 2°V[,g"d]e(e) + 2Q"j[e(e)Q^,]e(e)) , (66b) 

- /^(e) (°V, - r,) + \5Vg'f{^) Tf - tf{e) 0^(^)1 r,- (66c) 



U'Ae) 



Here R^^di^) is the curvature tensor of the metric (jabi^), and ^R'^bcd is the curvature 
tensor of the metric •jab- Note that ^R°'^cd = 7^^ ^R^cd is zero if at least one index is 
temporal, while if all indices are spatial 

^R\m = 2K6%6%. (67) 



^^See Bruni et al (1992) for a comprehensive treatment of linear perturbations using this formal- 
ism. 

^■^Malik and Wands (2004), equations (4.17) and (4.18) and Christopherson and Malik (2009), 
equations (4.11)-(4.13). 
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A. 2 First and second order gauge-variant perturbations 

Our first goal is to derive expressions for tlie geometric operators R"^cd (''''■'/) and 
T^°'\d{.fif) tliat determine tlie perturbed Riemann tensor tlirougli equations fl6al) . 
To accomplisli tliis we need tlie Taylor expansion of R°'\d{^)i given by ( l66l) . 

We begin by deriving expressions for ^^^g'^'^ and ^"^"fg"-^. We Taylor expand the 
relation 5% = g°''^{e) (jcbi^), which gives 



(Ikab _ _ fab {2)^ab _ _{2)eab _|_ 2 /ac rfe 


(68) 


where we used the following relations 




W(A5) = (0U«5 + (iU(%, 


(69a) 


(2\AS) = (°U (2)5 + 2 (lU (i)fi + (2U (°)5. 


(69b) 



We next Taylor expand (5V(e), as given by ( I64bl) . and use ( l69l) . This leads to 

^'~'Q%c = Q^(/), ('^QV = QV (('V) - 2/'^"Q,6c(/). (70a) 

where the operator Q%c is defined by 

Q^(/) := 7'^'(%/c)d - I ^dAc). (70b) 

Finally Taylor expanding W^^cdi^) as given by fl66|) and using fl6a|) . in conjunction 
with (1^ . (1^ and (I70ap . gives the leading order term 

R^U^^'V) = R^U^'V) + ^S^^iMi^^-^f), r = 1, 2, (71a) 

where 

R^U^^-^f) = -2 °V[e°V'" ('-^1 + (^^^/i'^ °i?''l'=.,, (71b) 

Q6^(Mj) ^ U^^,e(('V), (71c) 

and the source term 

n^'Uf. f) = ^"'cd(/, /) + ^6%l{%if, /), (72a) 

where 

TZ^Uf, f) = 2/i"(2R''l^,,(/) - //I °i?^4d) - 4Q^[%(/)Q|/P,](/), (72b) 

WV/,/) = -2//(u^V(/) + '^''/Q''^(/)^.)- (72c) 

In the above equationq^j we have also defined 

(^%{^''^f) := H'' (°Vd - r,) + 16", (^V^' r, + 7"^ Q%((''V)1 r/. (73) 



38To obtain ((7Tb)) we used ^Vf^OVrfj ('')/"'' = ^'''>fe^''°R''Kd, while OVc^'^V"'' = 2Q('^'')c(('^V) was 
used to obtain (j72bl) . 
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Equations ( 17T]) and ( 172]) constitute one of the main results of this paper. They 
express the first and second order Riemann tensor perturbations, as given by ( l6al) . 
in terms of the metric perturbations as given by equations ([3]), with the second order 
perturbation written as the sum of a leading order term and a source term. 

To proceed further we need to introduce local coordinates as in section 12. 2.1] 
which implies that the covariant derivative of a tensor, ^aA, and the gradient r^ 
Ia|) . assume the form 

'VoA = d„A, °V,A 



defined by 



D.A, 



ro = 71, 



0. 



(74) 



where d^j denotes partial differentiation with respect to r] and Dj is the spatial 
covariant derivative of ■jij. 

We have already used equations ( 17T1) to write the Einstein leading terms in the 
form ( IT9l) in the main text. We now use equations ( 1721) to write the Einstein source 
terms, as given by ( pT]) . in the form 



Gif, f) = -i [(1 + 34)7^ + 27^o + w\ - mcAUoJo)] 
g%{fj) = -2nB,ifoafon + Gl 



where 
>V.,(/,/) 

4(/,/) 

Q%fJ) 



An 

-2 



foaQ% — 2/ojDj/oo — /jfcQ j 



^"^(3(aj)kfa + ^'^^a{ifj)(3 " K{f\ifj)k 



]_ fk f 



V) 



'^Ja r< km + J^[J ml 



W'k 



TZi + 2K{^fQif fc — fomf"^i, 



\{f 



2Q"[\QrL, 



(/OO — 3/ fc)R ™am + /ofcR 



fcm 



m dOA; 



nR' 



4Q^[°„Q 



/3 m- 



(75a) 
(75b) 
(75c) 
(75d) 

(76a) 

(76b) 

(76c) 
(76d) 

(76e) 



For notational brevity we have dropped the arguments (/, /) of the quantities (1761) 
and the argument / of Q%c{f) and R"^cd(/) when they appear on the right side 
of (1751) and (1761) . It remains to give the components of R"''cd(/) and Q%c{f)^ that 
are defined by (I71bp and (]70bp : 

R%„ = \ (Dt + l^^mD^) /oo + 9,Q°^^, (77a) 

R%„ = 2D[,Q°^„], (77b) 

R%„ = -2D[*Ql°I^L + IKgSj, (77c) 

R*4„ = -2 (D[fcD[* + ir5[J*) /^f ] , (77d) 



Q° 



00 




-\d,f, 



00) 



Q Oi — ~oDi/o 



2^«J00) 
Q ij ~ 2^vfii ~ ^{ifj)0> QoO = f^r;/o* ~ 2^^ ^00, 

In summary, equations (175]) . in conjunction with (17^ -( 175]) . give the general expres- 
sions for the Einstein source terms. 



(78a) 
(78b) 
(78c) 
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A. 3 Gauge invariance and the Replacement Principle 

Gauge fields and gauge invariants 

In cosinological perturbation theory a second order gauge transformation can be 
represented in coordinates as follows: 

^--x'' + e^^^C + |e' (^'^r + ^'^C /'^^') , (79) 



X 



where '•^^^° and '•^^^^ are independent dimensionless background vector fields. Such 
a transformation induces a change in the first and second order perturbations of a 
tensor field A according to 

(i)i[e] = ('U + ^(i)^^")^, (80a) 

(2)i[e] = (^U + i:(.)^(°U + £a,^ (2(iU + i:(i)^(°U) , (80b) 

where £ is the Lie derivative (see, e.g., Bruni et al (1997), equations (1. !)-(!. 3)). 

One can impose restrictions on the tensor perturbations '■''•'A[,^] by letting ^'''\"' 
depend suitably on the dimensionless perturbations aJ'-^'^'^'A (r = 1,2), a procedure 
that can be referred to as perturbative gauge fixing. If the vector fields '^'^\"' up to 
order r are fully determined via the assumed perturbative restrictions, we say that 
the gauge is fully fixed to order r. An important special case of gauge fixing is 
order by order gauge fixing, in which the same conditions are imposed on the first 
and second (or higher) order perturbations. When the gauge has been fully fixed to 
order r, all remaining perturbative quantities to order r are rendered gauge-invariant. 

The key features of our dimensionless version of Nakamura's method for con- 
structing gauge invariants up to second order are as follows (see Nakamura (2007), 
and UWl, equations (5)-(8), for the linear case and further discussion and refer- 
ences). Given a family of tensor fields A{e) with a"'A{e) dimensionless we definqfj 

(i)A[X] := a" ((^U - £(i)x^°U) , (81a) 

(2)A[X] := a" pU - £c2)x^'^A - £,,,^ {2A^'^ - i^t.jx^^U)) . (81b) 

Comparing flHTj) and flHOl) multiplied by a", reveals that the equations have precisely 
the same form if we identify — ('^■)X" with ^'^\"-. Hence given a fully specified gauge 
choice ^'^\^, if we choose '^'^^X^ = —^'''\'^ then the ('')A[X,] will be gauge invariants 
that coincide with the a" '^'"•'yl[^,]. In other words, imposing conditions on the '^'")A[X] 
that fully determine the ('"^X'^ corresponds precisely to full gauge fixing. Due to the 
close relation between the roles of (''^X'^ and ^'^\"' we refer to (''^X'^ as gauge fields, 
and we say that '^-'^)A[X] and *^^-'A[X] are the first and second order dimensionless 
gauge invariants associated with ^^^A and ^^M, respectively, by X -compensation. 



^^Compare with equations (2.26)-(2.27) in Nakamura (2007) and (2.34)-(2.35) in Nakamura 
(2010). The factor of a" in our equations ensures that our expressions ('')A[X], r ~ 1,2, are 
dimensionless. 
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We note in passing that another way of ensuring that the ^^^A[ X] a re gauge- 
invariant expressions is to choose dimensionless fields ^"^^X"" that satisfjo 

^(i)X« = Wx^ _ (i)x« = (1)^", (82a) 

A(2)x" = (2)x'* - (2)x« = (2)^" + [(1)^, (i)x]". (82b) 

These conditions are obtained by applying flHOj) to f lHT]) and demanding that A(^)A[X] : 
and A(2)A[X] = 0. In UWl we made use of (I82ap in constructing gauge invariants 
(see section 2.2). At second order, however, the length of the computation makes it 
impractical to verify (]82bp directly. 

The one-to-one correspondence between full gauge fixing and the Nakamura ap- 
proach makes the choice between them a matter of aesthetics and personal prefer- 
ence o We choose to work in the formally 'dimensionless gauge-invariant '^''''A[X] 
picture', but everything we do has a 'fully gauge fixed picture' analogue (replace X 
with -^ and (''^aJx] with a" (^U[^]). 

In Appendix IA.4I we apply the Nakamura approach to the metric tensor to con- 
struct gauge fields and gauge invariants. 

The Replacement Principle 

The equivalence of gauge fixing and the Nakamura approach is expressed in the 
Replacement Principle. We present two versions, one for the Riemann and Einstein 
tensors, and one for the stress-energy tensor. 

The dependence of the perturbations of the Riemann tensor on the metric pertur- 
bations can be written symbolically in the form ( l6a|) . using the geometric operators: 

a^^'^R^^, = R-^^edl^'V), a'^'^R'^',, = R'^U'^'^f) + U^U'^'^f, '^'^f), (83a) 

where ^''^/ is shorthand for '•^■'/afe, with r = 1, 2. The Replacement Principle for the 
Riemann curvature states that the gauge invariants associated with '^'^'>R'^\d and with 
^^''fab by X- compensation are related by the same operators: 

(i>R'^^,4A] = K'^di!-^^), (%'^'',4X] = K'^di!-^ + 7^'^'ed(('^f , (')f), (83b) 

where *^''^f is shorthand for '^^^fabfA]. A similar result for the Einstein tensor can be 
derived from the above using ([5]). 

This Replacement Principle has its origins in the work of Nakamura, although he 
does not state it explicitly in the above form. See, for example, Nakamura (2010), 
equations (B9)-(B13), in conjunction with equations (2.36) and (2.37). 

The stress-energy tensor of a perfect fiuid, as given by fl50|) . can be viewed as a 
function of the variables F = {p,p,Ua,gab)- The perturbations of the stress-energy 
tensor can be written symbolically in the form: 



a' 



(^yr% = J%{^^^F), a^^^>r% = JM^^^F) + TM^^^F, (^^F), (84a) 



^'^These conditions arise in Nakamura's work. See for example Nakamura (2007), equations 
(2.23) and (2.25). 

^-'^The relation between gauge fixing and Nakamura's method has also been discussed by Christo- 
pherson et al (2011). 
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where T";, is the hnear leading order operator and T% is the quadratic source term 
operator, and ^'^^F = {^^^p, ^^)j, ^^\i'a, ^^^dab), with r = 1,2. The Replacement Principle 
for the stress-energy tensor states that the gauge invariants associated with ^^\r% 
and with ^'^^F by X-compensation are related by the same operators: 

^'>r%[x] = n((')F), (2>r"4x] = n((%) + r%i^^^, ^'^), (84b) 

where '•^^ is shorthand for ^^^[X]. This result can be deduced from Nakamura 
(2007)P 



A. 4 Construction of the gauge fields X°- 

We begin with the metric mode decomposition i^^. This decomposition assumes 
that the inverse operators (Green's functions) D~^, (D^ + 2K)~^ and (D^ + 3K)~^ 
exist (see UWl, and Nakamura (2007), page 19, for further discussion), as seen 
explicitly when one extract the various modes. In the present context this is accom- 
plished by means of the following mode extraction operators^ 

S' = D-^D*, S'^ = |D-2(D=^ + 3Ky^ B'\ (85a) 

V/ = (5/ - BiS^ , Vi^ = (D^ + 2K) "' Vi^D'^> , (85b) 

j:km ^ ^ik^rn) _ D(,V,f™ - D,,5'='". (85c) 

Applying these operators to (!24|) gives 

<^ = -|/oo, 5 = 57o^, B, = VUoj, (86a) 

C = -^S Jij, Cj = Vi Jjk, Cij = 2'ij J km, (8Db) 

V^ = -|(/-D25^V.,). (86c) 

Analogous relations hold for the modes of ^"^^fab- Explicitly extracting modes by 
means of the mode extraction operators fl85|) becomes essential at second order, and 
these operators are therefore used frequently in this paper. 

In order to construct gauge invariants associated with the metric perturbations 
we apply (ISTj) with ra = — 2 to Aab = gab = o,'^9ab and defincj 

iab[X] := fab - a~'^£{i)x{a^lab) , (87a) 

^X,[X] := ^%h - a-^£i2->j, {a^^ab) + J'abiX], (87b) 

where 

J'abiX] := -a~^£(i)x{2a^fab - <£(i)x («%«&))• (87c) 

We perform a mode decomposition of fafe[Xl and '^^^fa5[X] using equations (p^ as a 
model, and introduce an obvious notationQ This mode decomposition enables one 
to determine the gauge fields uniquely, as follows. 



"'^His equations (4.97)-(4.102) correspond to our equations (J84l) although it requires close scrutiny 
to conclude that his equations can be written in the operator form of (|84p . 

'*'^Note that D^ and D^- are defined in equation (I12|) . and that {ij) is defined in footnote fTTl 

^"'Here ^^^Aab — ^^^gab — a^lab, '''^''Aab — ^^hab — o? '^^^fab, and we denote the dimensionless gauge 
invariant associated with '^^^gab by ^''^fafc[X], r = 1,2. 

'^'^For example, ^^V ^ ^^^$[X], ^'^^B -^ ^'^'>B[X], as in equations 
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We begin by expressing fl87al) and (187bp in terms of the modes as defined by 
and the mode decomposed spatial vectors ^^^Xi = T)i^'^^X + ^^Xi, which gives 



^'^^[X] 

^''^C[X] 



0, 



(9, + ?^)Mx°-ij-oo[X], 



r)B + (0x0 - 5^(0x + S'J^oi[X], 
^H-9,Mx, + V/J-o,[X], 
0C<_Mx + l5'^'j-,^-[X], 

oc<^_Wx, + v/'=j-,4x], 

^'V + n^^'^x^ - liJ'Mx] - b^s'^:f,,[x]), 



(88a) 
(88b) 
(88c) 
(88d) 
(88e) 
(88f) 



where r = 1,2, and the source terms J-'ab[X] do not appear when r = 1. 

We next determine '^''-'X and ^^X^ uniquely by imposing the conditions (''''C[X] = 
and ('')Ci[X] = in dSSd]) and (l88e|) . respectively, which lead to 

(i)X = (ifc, (^)X, = (i)a, (89a) 

(2)x = (2)c + i 5^^- J-,,[X], (^)X, = (2)a + V/^j;-fc[X]. (89b) 

At this stage the mode decomposition for ^^^fab[X],r = 1, 2 assumes the forno 

Mfoo[X] = -2W<|.[X], (90a) 

Mfo^X] = D,MB[X] + (^)B,[X] , (90b) 

(^)f,,-[X] = -2M^[X]7,,- + 2(^)Q,[X] . (90c) 

It remains to determine ^''^X^, where (-"^^X" together with ^-^^X = ^^^C and ^^Xj = ^^^Ci 
form (^)X" which is to be inserted in ^ij [X] in (I89bp . There are two ways to determine 
^^^X^ algebraically in a unique manner. The first way is to set (^'^[X] = 0, which 
via f l88b|) with (El]) inserted yields 



«X0 = «Xp° :=-(% + 9, «C, 

(2)X0 = (2)XJ := -(2)5 + 9,(2)C - S''J^o^[Xp]. 



(91a) 
(91b) 



Note that (^^X^ together with fl89al) yields (^^X^. Substituting (^^X^J into P^j[X] and 
J-oi[-^p] in (!89bp and (I91bp . respectively, gives (^^Xp. We refer to the gauge fields 
(''•'Xp as the Poisson gauge fields since they result in the Poisson gauge invariants: 



(r)$ ._ (0$[Xp], (^)^ := ('■)^[X, 



pJ' 



^'^^^[X, 



PJ' 



(^)C,,[Xp], (92) 



when (1891) and (1911) are inserted into (IHHj) . 

The second way is to set '•^•*\E'[X] = 0, which via ( 88g ) with (I89ap inserted gives 



(i)X° = (1)X0 := -U-^^'^, 

(2)x0 = (2)x0 := _H-i(2V + |(n[X,] - D25''^-/-,,[XJ), 



(93a) 
(93b) 



*^This choice of the spatial gauge is essentially that made by Noh and Hwang (2004). See their 
equation (259) and the discussion on page 37. 
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where (^)X° together with fl89aD yields '■^^X^. Inserting '-'^^X^ into J'ij[X] in d89bl) 
and (l93b|) . and into J^\[Xc] in (IMbj) . gives (2)x^«. We refer to the gauge fields (^'^X;? 
as the uniform curvature gauge fields since they result in the uniform curvature 
gauge invariants^ 



WA:=W<l>[Xe], (^>B := (^>B[X,], (''^B^Xe], (^)Q,[X, 



(94) 



when ( !89l) and ( l93l) are inserted into ( l88l) . 

Note that at first order ^^^C^j[X] = ^^^dj, i.e., ^^^dj is independent of (^^X'*, 
while at second order (^)Cjj[X] depends on X, so that (^)Cjj[Xp] and (^)Cjj[Xc] are 
unequal. 

B Source terms in the uniform curvature gauge 

We make the replacement fab -^ ^abi^c] in the expressions f l75|) for the Einstein 
source terms using the metric perturbation ( 126|) in the uniform curvature gauge 
subject to the restriction ( ]33ap . i.e. 

foo[X,] = -2A, fo4Xj = D,B, f.,[Xj=0. (95) 

The constituent terms, as given by equations (176 1) -( ITHj) . can be evaluated separately. 
This calculation yields the following expressions for the source terms in the uniform 
curvature gauge, assuming a purely scalar linear perturbation: 

Qi.if,, fc) = 2 [2AD,, (A + CbB) + (9, A + iD^B) D,,B 

+ (D(,A)D,)(A + 2-HB) - (D^^.B) D'=,)B - ir(D(,B) D,)B] , 



(96b) 



Qii,, f,) = -1(1 + 3C^)7^ - lW\ + 2-H£^(-4A2 + (DB)^) 

- I [2AD2(A + CbB) + (9^A - 4HA) D^B + (DA)'] , 

^i(fe,fc) = D,7^-3?/^°, (96c) 

^° (f„ f,) = -2HDi(-4A2 + (DB)2) + g% (96d) 



where 



>V\ = An [2 A D'B + (D^^ A) D^B] , (96e) 

7^ = W\ - (D^„B) D™feB + |(D2B)2 - AK{BBY, (96f) 

g^i = -2(Dj A) D^iB + U'DiA) D^B - SKA D,B. (96g) 



We note that the source terms in the Poisson gauge as given by (l35l) were derived by 
imposing the restriction $ = \E', which led to significant simplification. If we impose 
the corresponding restriction in the uniform curvature case, namely A = —CbB (see 
equations (38a) and (40a) in UW2), then the above expressions simplify somewhat, 
as can be seen by inspection. 



^''We note in passing that that Christopherson et al (2011) have given the relation between the 
uniform curvature metric gauge invariants ^^^A/^^B and the Poisson gauge invariants (2)$ (2)v|/ (see 
their equations (4.56) and (4.58)). 
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